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Abstract 


This paper describes the design and evaluation of a stochastic opti- 
mal feed-forward and feedback technology (SOFFT) control architecture 
with emphasis on the feed-forward controller design. The SOFFT ap- 
proach allows the designer to independently design the feed-forward and 
feedback controllers to meet separate objectives and then integrate the 
two controllers. The feed-forward controller has been integrated with 
an existing high- angle- of- attack (high-a) feedback controller. The feed- 
forward controller includes a variable command model with parameters 
selected to satisfy level 1 flying qualities with a high-a adjustment to 
achieve desired agility guidelines , a nonlinear interpolation approach 
that scales entire matrices for approximation of the plant model, and 
equations for calculating feed-forward gains developed for perfect plant- 
model tracking. The SOFFT design was applied to a nonlinear batch 
simulation model of an F/A-18 aircraft modified for thrust vectoring. 
Simulation results show that agility guidelines are met and that the 
SOFFT controller filters undesired pilot-induced frequencies more effec- 
tively during a tracking task than a flight controller that has the same 
feedback control law but does not have the SOFFT feed- forward control. 


Introduction 

Typically, control designers try to achieve dif- 
ferent, and sometimes conflicting, design objectives. 
Examples of these control objectives include the fol- 
lowing: (1) good closed-loop stability characteristics 
(damping and bandwidth), (2) desired response to 
pilot input commands (flying qualities), (3) quick 
response to pilot input commands during large am- 
plitude maneuvers (transient response), (4) main- 
tenance of stability and tracking performance de- 
spite system uncertainties and various nonlinearities 
(robustness), (5) attenuation of high-frequency dis- 
turbances and measurement noise, and (6) accom- 
modation of low-frequency plant disturbances. 

Modern optimal control designs that attempt to 
meet many of these requirements with only a sin- 
gle cost function place conflicting demands on the 
controller. For example, when the single cost func- 
tion is optimized to provide good tracking charac- 
teristics, the controller usually has a high bandwidth 
with large feedback gains and poor noise attenuation. 

The SOFFT (stochastic optimal feed-forward and 
feedback technology) approach separates the feed- 
forward and feedback control objectives. The feed- 
back controller and the feed-forward controller are 
then designed with different cost functions. Control 
objectives such as 2, 3, and 4 that relate system re- 
sponse to pilot input commands are met by the feed- 
forward controller. Control objectives such as 1, 4, 
5 and 6 that relate to closed-loop damping and sta- 
bility, bandwidth, plant disturbance accommodation, 
and external noise reduction are met by the feedback 


controller. After the designs are completed, the feed- 
forward and feedback controllers are integrated with 
the SOFFT control structure to best meet all control 
objectives. A particular incremental implementation 
(described subsequently in this paper) is applied to 
the integrated controllers. This implementation al- 
lows control objective 6 to be met and avoids the 
need for operating point trim schedules. 

This paper describes the design for a SOFFT 
feed-forward controller that is integrated with an 
existing high angle-of-attack (high-a) feedback con- 
troller. The feed- forward controller includes (1) a 
variable command model with level 1 flying quali- 
ties and with a high-a adjustment to achieve agility 
guidelines, (2) a unique interpolation approach that 
scales entire matrices for approximation of the plant 
model, and (3) equations for calculating feed- forward 
gains developed for perfect tracking. The nonlinear 
batch simulation includes results for agility perfor- 
mance and tracking. 

Background 

Significant technical advancements in modern 
control theory have occurred during the past two 
decades. This section contains a summary of pre- 
vious research efforts that have influenced the direct 
digital control design described in this paper and the 
SOFFT (stochastic optimal feed- forward and feed- 
back technology) methodology (ref. 1). These previ- 
ous research efforts illustrate continuous technologi- 
cal improvements leading to a direct digital modern 
control design approach that is being used in a prac- 
tical airplane control design application. 



During the mid to late 1970’s, a digital flight con- 
trol system (refs. 2 and 3) was developed for use on 
the VALT (VTOL approach and landing technology) 
CH-47B research aircraft. The key technologies in- 
cluded a full-state direct digital design for inner-loop 
control, a PIF (proportional integral filter) control 
structure, and an incremental implementation for the 
digital flight computers. The main disadvantage of 
this technology was that the full-state control de- 
sign required state estimators for feedback. In addi- 
tion, complex gain-scheduling techniques were used 
to make the control system operational over a wide 
flight envelope. To demonstrate these methodolo- 
gies, a command generator and outer-loop controller 
were designed around the inner- loop PIF controller. 
The CH-47B aircraft then flew automatically along 
a curved four-dimensional trajectory to a hover and 
finally descended vertically to a landing (ref. 4). 

During the late 1970’s, the DIALS (digital inte- 
grated automatic landing system) control law was 
developed (ref. 5). DIALS was a full-state feedback 
design that required a full-state estimator. Because 
the application was for a limited flight envelope, only 
a single point design was used; thus, ad hoc gain- 
scheduling techniques that would have been required 
for a larger flight envelope were avoided. This con- 
trol methodology was flight tested on the Boeing 737 
research airplane at Langley during the early 1980’s 
(refs. 6 and 7). A constant-gain feedback matrix was 
used because automatic landing, which results in rel- 
atively constant flight conditions, was the overall con- 
trol objective of the flight test. Also, trim values for 
the steady-state operating condition were required 
because the control methodology was a full-state de- 
sign that was implemented in perturbational form. 

A discrete, optimal output, feedback algorithm 
was developed during the early 1980’s (refs. 8 and 9). 
This algorithm was an improvement over VALT and 
DIALS because the feedback controller only needed 
measured signals instead of the full-state feedback. 
A key advantage of this algorithm was that all 
important dynamics such as actuators, sensors, and 
filters could be included in the design process. A dis- 
advantage is that the optimal output feedback tech- 
nique was still a single point design approach, and 
ad hoc gain-scheduling techniques were still required. 
This technology, including the PIF controller with 
the incremental control structure, was used in re- 
structurable controls research applications (refs. 10 
to 12) in which control effectors were reconfigured 
to accommodate failures, and it was used in wind 
shear penetration research during thefcapproach and 
landing phases of flight (ref. 13). For all these appli- 


cations, nonlinear batch simulation of the Boeing 737 
at Langley was used as the test bed. 

During the mid 1980’s, an algorithm for multi- 
model output feedback was developed (ref. 14). The 
objective was to improve control design robustness by 
specifying different models around the same operat- 
ing point and by designing a fixed-gain output feed- 
back matrix that satisfied all models. With this al- 
gorithm, the control system should be robust enough 
to handle many different parameter variations. 

The multimodel methodology led to the next ma- 
jor advancement during the late 1980’s when the 
variable-gain output feedback methodology was de- 
veloped (refs. 15 and 16). This methodology allows 
simultaneous processing of multiple design points 
over the complete flight envelope, thus creating a 
more efficient design tool. The design algorithm gen- 
erates feedback gains that, in combination with a 
priori selected scalar design parameters, create an 
optimal gain schedule. The variable-gain technique 
is a modern control approach that can be used in 
a practical airplane control design application; how- 
ever, specifications for flying qualities cannot be di- 
rectly incorporated into the methodology. 

Variable-gain control was first applied to re- 
configurable aircraft (ref. 17), where gain- scheduling 
parameters were functions of the percentage loss of 
control effectors. Recent applications were for air- 
plane up-and-away flight (refs. 18 and 19). The 
feasibility of using the variable-gain output feed- 
back methodology was established by using four de- 
sign conditions and one gain-scheduling parameter 
for high-a flight at constant altitude (ref. 18). An 
expanded case with 39 design conditions and 6 gain- 
scheduling parameters was used to design a feed- 
back controller that covered the HARV (high-angle- 
of-attack research vehicle) flight envelope (refs. 19 
and 20). The HARV (ref. 21) is an F/A-18 aircraft 
that has been modified to include thrust vectoring. 
The controller described in reference 19 has been 
tested in real-time piloted simulation and is sched- 
uled for flight test. This controller is used in the 
SOFFT feedback system described in this paper. 

Some research on explicit model following is ger- 
mane to the SOFFT approach. Explicit model fol- 
lowing was one of the techniques investigated on 
an unstable advanced fighter model (ref. 22). A 
single quadratic cost function was used to mini- 
mize the error between the aircraft response and the 
command- model response. This controller had high 
gains and correspondingly high controller bandwidth, 
both of which resulted in excessive actuator rates. 
These results are typical of explicit model following 
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techniques that use a single cost function. An- 
other method used explicit model following but had 
separate designs for the feed-forward and feedback 
controllers (ref. 23). This approach allowed the feed- 
forward controller bandwidth to be tailored indepen- 
dent of the feedback controller bandwidth to meet 
flying qualities objectives. 

Finally, during the early 1990’s, the SOFFT tech- 
nology was developed (ref. 1). SOFFT is a design 
methodology that uses a unique control structure for 
integrating the feed-forward and feedback controllers 
and is applicable to multiple design points over the 
complete flight envelope. Also, flying qualities can 
be incorporated into the control design. A controller 
using the SOFFT methodology has been designed 
for application to a model of an HARV airplane. 
This controller incorporates many of the features de- 
scribed in previous references, such as (1) direct dig- 
ital design, (2) PIF control architecture, (3) incre- 
mental implementation, (4) variable-gain feedback 
design, and (5) separation of designs for the feed- 
forward and feedback controllers. This configuration 
was recently tested in a nonlinear batch simulation 
and is the subject of this paper. 
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Nomenclature n z 

Notation used in this paper includes boldface n ~ ,r 
symbols for matrices and vectors and italicized sym- n ZtSS 
bols for scalars. 


A; 

B z 

C n 


Cl 

C z 

c 

CAP 

ni 


e 


* 


command- model continuous-state 
matrix 

command-model continuous-control 
matrix 

normal-force coefficient 

variation of normal-force coefficient 
Cjv with a, rad -1 

interpolated plant-model state to 
output matrix 

command-model output matrix 

mean aerodynamic chord 

control anticipation parameter 

interpolated plant-model control to 
output matrix 

feed-forward command tracking 
error 


Ph 

P 

P 


Qc 

Q 

Q 

Q 

S 

s 

U X 


u 


* 

X 


H* plant-model state to regulated u z 

variable output transfer matrix 


matrix used to select regulated 
outputs 

height 

identity matrix 

command-model control input gain 

plant-model state gains 

integrator gain in feedback 
controller 

control filter gain in feedback 
controller, sec -1 

proportional feedback gain matrix 

command-model state gain matrix 

coefficient for sampling sequence at 
time 

number of plant models 

number of models used in scaling 
plant matrices 

number of gain-scheduling 
parameters 

normal acceleration, g units 

load- fact or command, g units 

steady-state normal acceleration, 
g units 

static pressure, lbf/ft 2 

parameter representing measured 
variables 

vector of measured parameters 

vector of parameters for each of the 
plant models 

impact pressure, lbf/ft 2 

pitch rate, deg/sec 

pitch acceleration, deg/sec 2 

dynamic pressure, lbf/ft 2 

reference wing area, ft 2 

Laplace variable 

control position command in feed- 
back controller, deg 

ideal plant-model control trajectory, 
deg 

command input to SOFFT feed- 
forward controller 
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V x 

W 

x* 

Vex 

yx 

yi 

Vx 

Vz 

z 

a 

a c 

r* 

r z 

AT 

A< 

Au z 

A x 2 

Ax* 

Ay* 

Ayi 

G 

0o 

p 

Pj 

<1>* 


rate command in feedback 
controller 

weight of airplane, lbf 

plant-model state vector 

command-model state vector 

ideal command to feedback 
controller 

output feedback measurement 
vector 

ideal plant output trajectory 

feedback controller error signal 

command- model output 

2 -transform variable 

angle of attack, deg 

angle-of- attack command, deg 

interpolated plant discrete control 
matrix 

command-model discrete control 
matrix 

discrete sampling period, sec 

incremental plant-model control, 
deg 

incremental cornmand-model input 

incremental command- model state 
vector 

incremental plant-model state 
vector 

incremental output feedback 
measurements 

incremental plant-model output 

pilot pitch stick input command, in. 

command-model short-period 
damping ratio 

initial pitch acceleration, rad/sec 2 

variable for gain-scheduling 
parameter selected by designer 

scaling parameter 

distance in parameter space from 
operating point to each plant model 

interpolated plant discrete state 
transition matrix 


command- model discrete state 
transition matrix 

interpolated matrix 

short-period frequency, rad/sec 

command-model short-period 
frequency, rad/sec 

Subscript: 

k coefficient for sampling sequence at 

time tfc 

Abbreviations: 

DIALS digital integrated automatic landing 

system 

HARV high-angle-of-attack research vehicle 

PI proportional plus integral structure 

PIF proportion integral filter 

SOFFT stochastic optimal feed-forward and 

feedback technology 

VALT VTOL approach and landing 

technology 

VTOL vertical take-off and landing 

SOFFT Feed- Forward Controller 

The complete feed- forward controller has two 
main components: the command generator with se- 
lectable modes and the SOFFT feed-forward con- 
troller (fig. 1). Pilot pitch stick input commands 6 sp 
go directly to the command generator, which scales 
the inputs based upon stick sensitivity, mode, and 
flight-operating conditions. Signals are calculated for 
two modes: a load-factor command mode n Z x and an 
angle-of-attack command mode a c . The mode selec- 
tor then chooses the smaller of the two signals for 
the output command u z ^. Selection of the small- 
est signal gives the best solution for the operating 
region of interest because each of the feed-forward 
gains can become large at different parts of the flight 
envelope. A derivation and implementation for the 
command generator with mode selection is presented 
in references 19 and 20. The methodology used in 
the references does not have a direct approach for in- 
corporating flying qualities guidelines. 

The SOFFT feed-forward controller converts 
command u z ^ into ideal trajectory commands for 
the control n* xk and ideal plant output y* fc , both 
of which are sent to the feedback controller. In 
this paper, the subscript k represents the coefficient 


<t>*x 

Wsp 
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Figure 1. Complete feed- forward control system. 



Figure 2. SOFFT feed-forward controller structure. 


for the sampling sequence at time and the aster- 
isk represents ideal trajectory variables with the as- 
sumption that the plant model represents the actual 
model. Figure 2 is a block diagram of the SOFFT 
feed- forward controller structure. The main com- 
ponents are the command model, plant model, and 
feed-forward gains. These components are described 
in the following subsections. 

Command Model 

The command model incorporates variable dy- 
namics, when* short-period frequency and damping 
are based upon the guidelines for level 1 flying qual- 
ities (ref. 24). An important flying qualities variable 
is the control anticipation parameter CAP , which is 
defined as the ratio of initial pitch acceleration 0\) to 
the steady-state normal acceleration n zss : 

CAP=A^^^]L (!) 

n z ,ss n z /n 

CAP is approximated as the ratio of the square of the 
short-period frequency uj sp to the ratio n z /a. The 
ratio riz/oc can be expressed as 

= C Nn qS = C Nn q 

a W 89.413 v 


where C\y o is the variation of the normal-force co- 
efficient C\ with a, q is the dynamic pressure, S is 
the reference wing area (400 ft 2 ), and W is the weight 
of the airplane (35 765 lbf). The calculation of CAP 
is used extensively at low o. but it is not a useful 
flying qualities metric at high a, where the airplane 
is in the stall region. Values of CAP that vary 
between 0.28 and 3.6 meet, the guidelines for level 1 
flying qualities. Figure 3 illustrates the relationship 
between CAP , u; iS p, arid the ratio n z ja and shows 
the region for level 1 flying qualities. The design 
parameters for the command model are located along 
the upper boundary of the level 1 region, with a lower 
boundary of 3 for the command- model short-period 
frequency lu z . The lower boundary of 3 was chosen to 
meet agility guidelines; related experimental results 
are shown later in the section entitled “Agility/’ 
Figure 4 shows a plot of C\ as a function of o 
with approximate values of slope C\r o located in the 
appropriate locations. The simulation model for C\ n 
is implemented by an upper boundary of 5.3 rad -1 
for a < 7.5°, a lower boundary of 2.4 rad' 1 for 
a > 17.5°, and a straight line between. The C\ 
was calculated from lift and drag coefficients, which 
were available at all trim conditions in our data base. 

At low q and high a (above 30°), lj z was modified 
to be 

w, = (0.6 + 0.08a) A?- (3) 

n z /a 
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Figure 3. Short-term pitch response command. (From ref. 24.) 



whereas the lower boundary can occur when either 
n z /a increases or a decreases. The smaller values for 
n z /a usually occur at high a where both Cjy tt and q 
are small. Equation (3) was derived experimentally 
to meet the high-a agility guidelines. Figure 5 is a 
plot illustrating this adjustment. 

Figure 6 shows the complete command-model 
configuration. Two items not previously discussed 
are the command-model short-period damping ratio 
( z and the discrete dynamic model. The value of ( z 
varies between 0.71 at high-speed conditions (n z /a 
= 10) and 1.0 at low-speed conditions ( n z /a = 1.1). 
The reason for this adjustment is because recent 
research shows that pilots prefer high damping at 
high-a flight conditions (ref. 25). 

The command model has second-order dynamics 
as shown by the transfer function 

Vzk. _ ^3 (4) 

u zk s 2 + 2 C, z u z +ujI 

where the variables oj 2 and Cz vary with parameter p 
(representing a and q in this case) as flight conditions 
change. Equation (4) is transformed to a discrete 
state-space representation 

x z (p,k + 1) = & z (p)x. z {p,k) + r z (p)u z (p,k) ) 

\ (5) 

Vz{p, k) = C z(p)xz(p, k) ) 


5 i— 


o 4 
<u 

C/5 


3 3 


0 30 


35 


40 45 

a, deg 


50 


n~ 



55 


60 


where x z is the command-model state vector, u z is 
the command- model input (single command), 4> z is 
the command-model discrete state transition matrix, 
T z is the command- model discrete control matrix, C z 
is the command-model output matrix, and integer k 
represents the sample at time t^. Matrices <& z (p), 
and r z (p) must be calculated at each iteration ( C z 
is a constant matrix for the selected implementation) 
and are approximated as 

&z(p) = exp[AT A z (p)} « 1 + AT A z {p) 

(AT) 2 A 2 (p) 


+ 


(6) 


Figure 5. Low-speed high-a adjustment. 

with a lower boundary of 3 and an upper boundary 
of 5. As an example, the upper boundary can occur 
when a is 55° or greater and n z /a is 1.2 or less, 


r z ( P ) = Ap(p) [<Mp)-i]eMp) (7) 

where A z (p) and B 2 (p) are the continuous matrices 
corresponding to the transfer function in equation (4) 
and AT is the discrete sampling period (0.0125 sec). 
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Because AT is small, a second-order approximation in equation (6) gives sufficient accuracy. Individual elements 
for equations (6) and (7) are 


$ = 


, (A Tw *) 2 

1 2 


AT(1 - ATQ z lj z ) 


L ATcuf(-l + ATQzUiz) 1 - 2AT CzOJz + 2(AT Qu) z ) 2 - (A7 ^ 


T, = AT Jl 


0.5AT 
1 - ATuzQ 


( 8 ) 

( 9 ) 


where these elements are calculated at each time iteration and argument p is neglected for simplicity. 



Plant Model 

The plant model produces the trajectory com- 
mands y* k and u* A . for input to the feedback con- 
troller. To produce reasonably accurate commands, 
the plant model must be approximated at any flight- 
operating condition. Figure 7 shows the plant-model 
configuration and illustrates how flight-measured pa- 
rameters are used to obtain matrices for the plant- 
model dynamics. This section described the plant- 
model interpolation and the plant-model dynamics. 

Plant-model interpolation . In general, plant 
models can be composed of a series of models repre- 
senting components such as actuator dynamics, air- 
plane dynamics, sensor dynamics, and filter dynam- 
ics. Airplane dynamics change continuously with 
flight-operating conditions. To accommodate the 
complete operating range, interpolation between a fi- 
nite number of specified plant models is necessary. A 
unique interpolation approach that scales complete 



Figure 7. Plant -model configuration (perfect sensors 
assumed). 

matrices is presented in reference 1 and is presented 
here for completeness. Linear interpolation is also 
a feasible approach, but it was not used because 
of time constraints and the additional complication 
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when more than two independent variables are used. 
Precise accuracy is not required for the feed-forward 
controller. Thus, compared with linear interpolation, 
the method used here is a relatively fast procedure. 

A finite number of plant models are specified as 
a function of several parameters p that vary over 
the flight envelope. A metric pj, representing the 
distance in parameter space from the operating point 
to each design model location, is computed as 

Pj(p) = ||p-P j || 2 = (10) 

where p is a vector of measured parameters that 
can be either linear or nonlinear, and p- 7 is a vector 
of equivalent parameters for each of the M plant 
models. All pj are sorted from smallest distance 
to largest distance and then the smallest n numbers 
arc selected to be used. A scaling parameter p(p) is 
calculated by using the n closest models as follows: 


pip) 


i 


n 


E 

i - 1 


I 

Pi(p) 


( 11 ) 


state vector; y* is the ideal plant-model output 
vector; and u* is the ideal plant-model control (scalar 
for the controller in this paper). For implementation, 
equations (13) are solved in incremental form where 
the input is Au*, which is defined as the difference 
between the values at two successive sampling times: 

A u*(p, A) = u* (p, k) - u* (p, k - 1) (14) 

and the output is A y*. Using the incremental im- 
plementation eliminates trimming problems because 
the incremental plant-model states are always zero 
during steady-state conditions. For the incremen- 
tal implementation, the dynamics are assumed to be 
constant during each sample interval. The total out- 
put is then solved by accumulating all previous in- 
crements as follows: 

y *(p, k) = y*.(p, A’ - 1) + A r/*(p, k) (15) 


Feed-Forward Gains 

Reference 1 presents an optimal cost function that 
is quadratic in states, controls, and the feed-forward 
tracking error which is defined as 


and the ratio p{p)/pi(p) is then used to weigh the 
various matrices of the n closest models as follows: 

«»> = £ ^ ( 12 ) 

where 0* (p) represents the interpolated matrix to be 
used for the plant model. If pj(p) equals 0 for any 
design model, then that model is used as the plant 
model and equations (11) and (12) arc not needed. 

In figure 7, p is a function of a, impact pressure 
Qc< and static pressure P s , in the design example for 
plant-model interpolation where each p* is limited to 
values within the design envelope. The design exam- 
ple used 39 plant models (M = 39) and 3 models 
(n = 3) for the interpolation process. The selection 
of n = 3 is a judgment factor based upon some pre- 
liminary analysis of the interpolation error. 

Plant-model dynamics . The plant model is 
solved as a discrete state-space representation as 
follows: 


e*{p,k) = Hy y*r(p, k) -y z {p,k) (16) 

where H ?y is a matrix that allows a selected combina- 
tion of plant-model outputs to track the command- 
model output at every instant of time. Two ap- 
proaches for gain calculation are shown in reference 1. 
The first is a variable-gain approach (refs. 15, 16, 
and 18), and the second is a perfect tracking ap- 
proach. The variable-gain approach generates an op- 
timal gain schedule in which the gains are optimized 
over the entire flight envelope and are functions of 
measured parameters at each instant of time. In the 
variable-gain approach, there is a trade-off between 
feed-forward state variations, control variations, and 
tracking performance. 

In this paper, the perfect tracking approach is 
used; that is, the optimal cost-function penalty 
weights on states and controls are zero. Feed-forward 
gains are generated to make c£ zero at all instants of 
time. The control law is in the following incremental 
form: 

A u*{p, k) = -K{p) Ax* (p, A) - K z (p) A x*(p, A) 


x* (p, k+\) = &l ( p ) x* (p, A:) + r* (p) u* (p, k) 1 

\ ( 13 ) 

y*(PiA) - C*(p)x*(p, k) + D*(p) u*(p, A*) J 

where 3>*, T*, C*, and D* are interpolated matrices 
that are updated each iteration; x* is the ideal plant 


- K u (p) Au z (p, k) (17) 

where K* is the plant-model state gains, K z is 
the command- model state gains, and K u is the 
command-model control input gain. (See fig. 2 for 
structure.) Equations (5) and (13) with D* assumed 
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zero, equation (16) with e£ equal to zero, and equa- 
tion (17) give the perfect tracking feed- forward gains 
as 


K n (p) = -[H*.(p)T*(p)}- 

] C Z (I>)T : (P) 

(18) 

k : ( p ) = -\r*(p)K(pT 

'c M*z{p) 

(19) 

K*(p) = -[K(p)K(p)]~ 

' H» 

(20) 

with 



H» = H„ C» = [1 

i i]c» 

(21) 

Equation (21) shows that the plant-model output 
is the sum of three regulated outputs that have 
been chosen to match the feedback controller. The 
variables in equations (18) to (21) are from the plant 


model and command model with the feed- forward 
gains calculated at each sampling interval. When 
D* is not zero, as with acceleration output., there 
is a small residual tracking error. However, D* 
is included in the dynamic equations for the plant 
model (eq. (13)). Because there is one feed-forward 
control, the matrix inversion in equations (18) to (20) 
is trivial. For two or three controls, the matrix 
inversion can still be accomplished reasonably fast. 

Algorithm 

The order in which equations are being imple- 
mented has been changed slightly from that shown 
in reference 1. The three main changes are as 
follows: 

1. The command-model input uses total variables 
instead of incremental variables; thus, the output 
has total variables. 

2. The total command-model output variable in- 
stead of the plant-model output variable is passed 
to the feedback controller. 

3. The incremental plant-model control signal in- 
stead of the total plant-model control signal is 
passed to the feedback controller. 

The first change was made because at the com- 
mand model output during the accumulation process, 
the incremental version resulted in a significant error 
caused by the variable dynamics, which are assumed 
to be constant during each sampling interval. The 
second change was made to eliminate the need for 
an accumulator at the plant-model output because 
perfect tracking is used. An accumulator acts as 
a summing device and is similar to an integrator. 
This approach also eliminates the small error caused 
by the nonlinear implementation because D* in the 


plant model is nonzero. Finally, the third change 
incorporates a limited accumulator in the feedback 
controller: thus, the need for a separate limited accu- 
mulator in the feed-forward controller is eliminated. 

Starting with input u z (p, k), the algorithm is as 
follows: 

1. Solve for the incremental input Au s (p, k) = 
u z (p. k) — u z (j) , k — 1 ). 

2. Compute & z (p) in equation (8) and T z (p) in 
equation (9). (Note C z (p) is a constant matrix 
in this paper.) 

3. Update the command-model dynamic equations 
(eqs. (5)). 

4. Compute the incremental command- model state 
vectors A x c (p, k) — x z (p. k) - x~(p, k — 1). 

5. Perform plant-model interpolation bv using equa- 
tions (10) to (12) for <£*(/)). r*(/>). C*(p), and 

d;(p). 

6. Compute feed-forward gains with equations (18) 
to (21 ). 

7. Solve' for the incremental control (eq. (17)). 

8. Update the plant-model dynamic equations 
(eqs. (13)) by using the incremental control in- 
put from step 7. (Note, the states and outputs 
are also incremental.) 

9. Perform either step a or step b. 

a. Solve equation (15) for y x (p<k) and compute 

y*(p< k) as in equation (16). 

b. Use y z (p<k) for H y y*{p.k) since for perfect 
tracking c*(p, k) is 0 (eq. (16)). The signal 
V*r(P'k) can represent either the command- 
model output or the plant-model output. 
(This step was used in this example.) 

10. Send variables A u*(p, k ), y* x (p. k)< and Ay *.(;>, k) 
to the feedback controller. 

Feedback Controller 

Design of the feedback controller is described in 
references 18 and 19 and is also summarized in this 
section for completeness. Figure 8 shows the discrete 
PIF control structure that is used for design and 
linear analysis. The output feedback measurement 
vector y xik , is input to both the proportional plus 
integral (PI) feedback paths, which are in parallel. 
The outputs from the PI structure go to a low-pass 
filter to produce the control position output u r ^. The 
PIF controller is a rate-command structure where 
the proportional feedback gain matrix J£ xx (p)* the 
integrator gain K x f(p). and the control filter gain 
K xu {p) all join at a summing junction to produce the 
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Figure 9. Feedback controller implementation. 


rate command v x k- The time step AT from v x ^ to 
u x k accommodates necessary computational delays. 
The regulated variables are defined by row vector H y 
as given in equation (21). 

The feedback gain matrix has a linear, func- 
tional relationship with parameter p and contains 
both constant- and variable-gain parts that are im- 
plemented as 

rip 

K( P ) = K 0 + X>(^) K.; (22) 

1 = 1 

where the variable represents some measured vari- 
able that the designer selects for the gain-scheduling 
parameter and n p is the number of gain-scheduling 
parameters. The relationship between p % and can 
be either linear or nonlinear. The feedback gains in 
figure 8 are partitions of the overall gain matrix as 

K(p) = {K xx (p)K xI (p)K xu (p)} (23) 


The feedback controller was implemented incre- 
mentally with total measurable quantities. (See 
refs. 1, 3, and 26.) The advantage of the incremental 
implementation is that trim tables are not required; 
thus, the airplane automatically goes to a new equi- 
librium state as the integrated output follows the 
command. Because this controller is described in ref- 
erence 19, only a few key equations are included here 
to show how the feed- forward signals are integrated 
into the feedback controller. Integration of the two 
controllers was found to be easy. (Some of the no- 
tation used in this paper differs from that used in 
ref. 19.) 

As shown in figure 9, the rate command v x {p,k) 
is solved as 

fix (p, k) = [I - AT K xu {p)\ v x (p, k - 1) 

- K jjr (p) [Ajte(p, k) - A y*{p, fc)] 

- A TK xI (p) [H y y x (p,k- 1) - y‘„(p,k - 1)] (24) 
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Figure 10. Major dimensions of HARV. All dimensions are in feet. 


Then, u x (p.k) is calculated as 

Mx(Prk) - u x (p,k - 1) 4- A u* r {p,h) + A Tv x (p,k - 1) (25) 

where A y x (p,k) is the incremental output feed- 
back measurement vector. Equation (25) shows 
that u x (p,k) is a combination of the integral of 
v x (p,k) (delayed one time period) and the accumu- 
lated sum of the feed-forward incremental control sig- 
nal Au*(p,k). 

Facilities 

The SOFFT controller was applied to a nonlinear 
simulation model representing an F/A-18 airplane 
modified to have multiaxis thrust vectoring for ad- 
ditional pitch and yaw control. This modified con- 
figuration is known as the HARV (ref. 21). The 
F/A-18 airplane is a multirole fighter/attack airplane 
with supersonic dash capability and good low-speed 
high-a maneuvering capability. Major dimensions of 
the HARV are shown in figure 10. Thrust-vectoring 
capability was added to the basic F/A-18 aircraft 
by removing the secondary nozzles and adding three 
thrust- vectoring vanes per engine (fig. 11). The 
modified airplane has a nominal weight of almost 
36000 lb, which is approximately 4000 lb heavier 
than the basic F/A-18 aircraft. 

The F/A-18 aircraft simulation on which the 
HARV model is based is discussed in detail in ref- 
erence 27. The HARV simulation was built from 
fully nonlinear aerodynamic, engine, and control sys- 


tem models of the production F/A-18 aircraft; these 
models were obtained from McDonnell Douglas Cor- 
poration. The McDonnell Douglas aerodynamic data 
base is for a = -10° to 90°, sideslip (3 = -20° to 
20°, altitudes to 60 000 ft, and speeds to Mach 2.0. 
Aerodynamic increments were added to the database 
because of the addition of thrust-vectoring vanes, ac- 
tuator housings, and a spin parachute. Jet-induced 
effects were added for the change in airflow over the 
airframe that resulted from thrust vectoring. 

The engine model, also obtained from McDonnell 
Douglas, incorporated thrust vectoring, the effects 
of Mach and altitude, and the dynamic response of 
engine thrust. Also included were the effects of a 
and vane deflection. Gross thrust and ram drag were 
tabulated separately; this tabulation allowed thrust 
vectoring to # act on gross thrust only. 

The SOFFT longitudinal control law discussed 
herein was integrated with a high-a lateral- 
directional controller to provide stability and maneu- 
verability in the lateral-directional axes. The thrust- 
vectored outputs from the two controllers go to a 
vane control system know as the mixer/predictor. 
The mixer/predictor converts pitch, yaw, and roll 
thrust vectoring commands into equivalent com- 
mands for the six thrust-vectoring vanes to yield the 
required jet deflection. 

Simulation Results and Discussion 

The SOFFT feed-forward controller was inte- 
grated with the HARV longitudinal controller 




Figure 11. HARV with thrust vectoring vanes. 
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(ref. 19) by inserting SOFFT between the command 
generator and feedback controller (fig. 1) and by elim- 
inating the command generator tracker. This config- 
uration was then tested in a nonlinear batch simu- 
lation to evaluate agility performance and tracking. 
Some specifics related to the feedback controller de- 
sign and the results of the nonlinear batch simula- 
tions are described in this section. 

Feedback Controller Specifics 

In the variable-gain feedback design, six gain- 
schedule parameters p r ( Vj ) were used. These param- 
eters are functions of o, Q ( and P s and were selected 
to cover independent degrees of freedom (a, speed, 
and altitude). The Pi(i/j) and their limits were 


pi = 0.1a 

(1.5 < a < 65) 1 

p2 = O.OlQe 

(10 < Q c < 470) 

P3 = 0.001 P, 

(498 < P s < 1200) 

P4 - Qc/Ps 

(0.008 < p 4 < 0.4) 

p 5 = 0.1a -3.5 

(a > 35) 

P 5=0 

(a < 35) 

P0 = 0.01Q r - 2.5 

{Qc > 250) 

P6 = 0 

{Qc < 250) ; 


where the limits were selected to cover the HARV 
flight envelope. These p, are used in equation (22) to 
calculate new feedback gains at each sample time. 

The feedback gains changed continuously with 
the measured variables. The function was completely 
continuous and smooth except at two points. The 
first four parameters cover the entire HARV flight 
envelope; the last two parameters cover only portions 
of the envelope and were selected to tune in those 
portions of the design envelope. Parameter p$ was 
used only when a is 35° or greater, and parameter 
Pq was used when Q ( is 250 lb/ft 2 or greater. Both 
P 5 and p$ have lower limits of zero and are not 
differentiable at the break points. When any value of 
vi exceeded the design limit, the variable was limited 
at the value shown in equation (26). The coefficients 
were selected to keep the Pi(uj) near unity. 

Although it is not shown in figure 9, u x {p, k) 
splits into two commands. One command is position 
limited and goes directly to the stabilator input, 
while the other command passes through a limited 
washout filter and becomes the pitch thrust- vectoring 
command. 

Agility 

The simulation approach for agility performance 
evaluation is to first trim the airplane at a pre- 
determined a (or Mach number) at 1 g and 25 000 ft. 



Figure 12. Pitch rate agility for various initial o trims at 
25 000 ft. 



ex. deg 

Figure 13. Pitch acceleration agility for various initial a trims 

at 25 000 ft. 

At low-a conditions, the thrust available was suffi- 
cient to maintain a flight path angle of 0°. How- 
ever, at a flight path angle of 35° or greater, max- 
imum thrust was insufficient and resulted in a neg- 
ative flight path angle. At time 0 + , the throttle is 
moved to maximum (afterburner), 2 sec later full aft 
pitch stick (5-in. maximum) is applied, and the peak 
pitch acceleration and pitch rate are measured. 

Figure 12 shows the pitch rate agility for various 
initial trim a’s, and figure 13 shows the pitch acceler- 
ation agility for the same conditions. In each figure, 
the solid line represents the guideline for desirable 
agility response (ref. 28) and the dashed lines rep- 
resent controller performance. The two data plots 
deviate above a = 30° where the high- a adjustment 
discussed previously is the only difference between 
the data of the command model with lower limit of 3 
for ujy and no adjustment and the data for the com- 
mand model with the low-speed high-a adjustment 
for uj z . The adjustment improves both pitch rate, 
q , and pitch acceleration, q , agility. Pitch rate ei- 


ther exceeds or meets the guideline at all a up to 
45°, when the data becomes slightly lower than the 
guideline. Pitch acceleration with the high-a adjust- 
ment exceeds the guideline at all a and significantly 
exceeds the guideline at most a. 

Tracking 

Approximately 19 see of pilot pitch stick input 
command 6 sp were extracted from a real-time simula- 
tion session during a tracking task for a different con- 
troller at approximately Mach 0.4 and an altitude of 
25 000 ft. These data were then used in the nonlinear 
batch simulation to evaluate the SOFFT controller. 
Shortly after the start of the simulation run, the pilot 
slowly increased throttle to afterburner and reached 
maximum thrust in approximately 4 see. The pilot 
also moved the lateral stick to roll the airplane to 
approximately 60° within the first 5 sec, and then 
maneuvered between 60° and 80° for the remainder 
of the simulation. Because this paper is for a lon- 
gitudinal controller, only those variables relating to 
the longitudinal axis are described. 

Figure 14 shows 10 time histories: 6 sp , u^, y* x ^ 
u) z , Gi c£, a, q , y T k' and u x f. where arguments p and 
k are dropped for simplicity. The last four time histo- 
ries (a, q , y :r j., and u x i) show a comparison between 
the SOFFT controller response and the controller re- 
sponses from reference 19. The first two time histo- 
ries {t) sp , u z k) are identical for both controllers, and 
the third to sixth time histories (z/*^.* G? c£) only 
apply to the SOFFT controller. 

Signal u z k has the same noise content as 6 sp , but 
it is larger in magnitude because of the command 
generator gains. Signal y* r ^ (step 9b of algorithm) 
illustrates the filtering within the command model. 
The plot of y* xP is clearly less noisy than that of 
particularly at the high-a flight conditions, where u) z 
is at or near the lower boundary of 3, and G ls near its 
upper boundary of 1. The feed-forward tracking error 
c p is relatively small (generally less than 0.1 percent 
of 1&*) and is caused by the assumption that D* 
(eqs. (13)) is 0 in the perfect tracking equations. 

Because the feedback controller regulates the sum 
of measured signals a, q and n z , a only approximates 
Perk' ^ 1C nonlinear simulation, pitch rate output 
is modified (not shown in fig. 9 for simplicity) by 
nonlinear gravity compensation prior to use by the 
feedback controller. The feedback controller tracking 
error y x ^ is defined as 

Vxk = H '/ y.ii - y* T .k = (o + q + n z ) k - y* Tk (27) 

which is the integral term (one time step ahead) 
in equation (24). As shown in figure 14, y x has 
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Figure 14. Tracking simulation. Mach = 0.45: h — 25 000 ft. 


peaks of 1.6 for SOFFT and is generally below 1.0, 
which indicates improved feedback regulation com- 
pared to the reference 19 controller which has peaks 
that reach 7.5. Comparison between the SOFFT 
controller and the controller in reference 19 shows 
that SOFFT filters undesired pilot-induced frequen- 
cies much more effectively, which results in signifi- 
cantly reduced control-effector activity u x j-. 

Concluding Remarks 

This paper describes the design and evaluation 
of a stochastic optimal feed-forward and feedback 
technology (SOFFT) control structure with empha- 
sis on the feed- forward controller. The feed-forward 
controller was designed separately with an objec- 
tive of perfect tracking, and then it was easily 
integrated with a previously designed feedback con- 
troller with different objectives. The main compo- 


nents of the feed-forward controller are the com- 
mand model, plant model with interpolator, and 
feed- forward gains. 

The command model incorporates variable dy- 
namics for the command- model short- period fre- 
quency ( LO z ) and the command-model short-period 
damping ratio (( 2 ), which are functions of angle of 
attack (a) and dynamic pressure ( q ). In general, as 
speed increases at low a, io z increases and £ z goes to 
its lowest limit of 0.71. Design parameters have been 
chosen along the upper boundary of the level 1 flying 
qualities guidelines with a minimum uj z of 3 to meet 
agility guidelines. At low dynamic pressure and high 
a (a > 30°), an adjustment was made to vary u) z as a 
function of a and q , with an upper limit of 5 to meet 
the high-a agility requirements. The largest damp- 
ing ratio of 1 occurs during low-speed flight, which 
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is generally at high a because pilots prefer greater 
damping during this phase of flight. 

A unique interpolation is used within the feed- 
forward controller to generate a plant model from a 
finite number of state-space design models that cover 
the flight envelope. A metric representing the dis- 
tance in parameter space from the operating point 
to each design model is computed and used to scale 
an entire matrix. Precise accuracy is not required 
for the feed-forward task, so the main benefit is 
that this procedure is relatively fast compared with 
linear interpolation using three independent vari- 
ables. Thirty-nine reference models were used in 
the design, and the three reference models closest to 
the measured operating condition were used for the 
interpolation process. 

A perfect tracking algorithm was used for on-line 
calculation of the feed-forw r ard gains. This algorithm 
was derived by solving for the optimal gains that 
make a selected combination of the plant-model out- 
puts follow the command model output. The equa- 
tions should be easy to calculate for as many as three 
controls. Integration of the feed-forward controller 
with the feedback controller by using the SOFFT 
structure and an incremental implementation was 
straightforward. 

Nonlinear batch simulation results show that use 
of the SOFFT controller enables agility guidelines 
for pitch rate and acceleration to be met. Without 
the high-Q adjustment, pitch rate agility was slightly 
below the guideline for a > 30 °. 

Tracking task time history plots comparing the 
SOFFT controller with another controller with the 
same feedback system shows that SOFFT filters un- 
desired pilot input frequencies much more effectively, 
has a smaller tracking error, and has reduced control- 
surface activity. Real-time and full-scale flight test 
pilot evaluations are still needed to determine flying 
qualities and tracking performance. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
March 31, 1994 
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the plant model, and equations for calculating feed-forward gains developed for perfect plant-model packing. 
The SOFFT design was applied to a nonlinear batch simulation model of an F/A-18 aircraft modified for 
thrust vectoring. Simulation results show that agility guidelines are met and that the SOFFT controller filters 
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